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NLO QCD corrections to pp — > ttbb + X at the LHC 
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We report on the calculation of the full next-to-leading order QCD corrections to the production 
of ttbb final states at the LHC, which deliver a serious background contribution to the production 
of a Higgs boson (decaying into a bb pair) in association with a tt pair. While the corrections 
significantly reduce the unphysical scale dependence of the leading-order cross section, our results 
predict an enhancement of the ttbb production cross section by a if-factor of about 1.8. 

PACS numbers: 14.65.Ha 



Extending earlier work where we discussed the 
next-to-leading order (NLO) QCD corrections to ttbb 
production via quark-antiquark annihilation, in this let- 
ter we present first results on the full NLO QCD correc- 
tions to pp — > ttbb + X at the LHC, i.e. we complete the 
existing results by the contributions from gluonic initial 
states. 

The QCD-initiated production of ttbb final states rep- 
resents a very important background to ttH production 
where the Higgs boson decays into a bb pair. While early 
studies of ttH production at ATLAS and CMS 0] sug- 
gested even discovery potential of this process for a light 
Higgs boson, more recent analyses 0, 0, @, 0] with more 
realistic background assessments show that the signal sig- 
nificance is jeopardized if the background from ttbb and 
tt + jets final states is not controlled very well. The cal- 
culation presented in this letter renders improved signal 
and background studies possible that are based on NLO 
predictions for ttbb final states. NLO QCD corrections 
are already available for the ttH_ signal jg, |jj fie] . [ T ill and 





FIG. 1: Some generic hexagon diagrams for hadronic ttbb 
production at NLO QCD. 
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the background from tt + jet 

On the theoretical side, the calculation of NLO cor- 
rections to 2 — > 4 particle processes represents the cur- 
rent technical frontier. The complexity of such cal- 
culations triggered the creation of prioritized experi- 
menters' wishlists lH, HI for missing NLO calculations 
for LHC physics, and the process of ttbb production 
ranges among the most wanted candidates. In recent 
years the field of NLO calculations to multi-particle pro- 
cesses received an enormous boost, most notably by ad- 
vanced methods for evalua ting one-loop tensor integrals 
for Feynman diagrams [H ESE IM III HI IH and by 
new methods employing unitarity cuts of one-loop am- 
plitudes analytically (see, e. g., Ref. 12411 and references 
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therein) or numerically [25 




numerical unitarity-based approaches have successfully 
passed their proof of principle in the calculation of spe- 
cific one-loop QCD amplitudes, including in particular 

multi-gluon amplitudes m eh, Hi], ud -> w+^g [m, 



W_ + g_gg [33. l34j. uu/gg -> ttbb, uu -> W+W~bb, 
uu — * bbbb, uu/gg — ► ttgg and in the evaluation 
of the leading-colour contribution at NLO to W+3jet 
production at the Tevatron P, [H, El. The Feynman- 
diagrammatic approach has already been used in com- 
plete NLO predictions for some 2^4 reactions at 
e + e~ 37|, [38| and 77 colliders [39l an d in the evalua- 
tion of the amplitude of qq — ► bbbb [40| at one loop. The 
evaluation of ttbb production — also based on Feynman 
diagrams, as documented in Ref. [l[ and this letter — 
represents the first full NLO calculation for a 2 — > 4 
process at a hadron collider. 

In LO QCD, 7 and 36 different Feynman diagrams 
contribute to the production of ttbb final states via qq 
annihilation and gg fusion, respectively. The virtual 
QCD corrections comprise about 200 one-loop diagrams 
for the qq and about 1000 diagrams for the gg initial 
state, the most complicated being the 8 and 40 hexagons 
for the respective channels. Some hexagon graphs are 
depicted in Figure [TJ The real QCD corrections com- 
prise gluon bremsstrahlung in the qq and gg channels, 
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qq/gg — > ttbbg, and (anti)quark-gluon scattering pro- 
cesses, Vg — * ttbbg.' We consistently neglect contri- 
butions involving b quarks in the initial state because 
of their suppression in the parton distribution functions 
(PDFs). In the following we briefly describe the calcula- 
tion of the virtual and real corrections. Each of these con- 
tributions has been worked out twice and independently, 
resulting in two completely independent computer codes. 

The evaluation of the virtual corrections starts with 
the generation of the one-loop amplitudes via two in- 
dependent versions of FeynArts jU, Usin g ei - 
ther inhouse Mathematica routines or FormCalc j43|. 
each diagram is decomposed in terms of standard spin 
and colour structures, as described in Ref. [1] for the 
qq channel in detail. The coefficients in the resulting 
linear combination of these standard structures contain 
the one- loop tensor integrals. The obtained expressions 
are not reduced to standard scalar integrals analytically, 
but the tensor integrals are evaluated by means of al- 
gorithms that perform a recursive reduction to master 
integrals in numerical form. This avoids a further in- 
crease of the huge analytic expressions and permits to 
adapt the reduction strategy to the specific numerical 
problems that appear in different phase-space regions. In 
detail, the 6-/5-point integrals are directly expressed in 
terms of 5- /4-point integrals ID, 2^] ■ Tensor 4-point and 
3-point integrals are reduced to standard scalar integrals 
with the Passarino-Veltman algorithm [3] as long as no 
small Gram determinant appears in the reduction. If 
small Gram determinants occur, the alternative schemes 
of Ref. Q are applied. Ultraviolet (UV) divergences 
are regularized dimensionally throughout, but infrared 
(IR) divergences are treated in different variants, which 
comprise pure dimensional regularization with strictly 
massless light quarks (including b quarks) and a hy- 
brid scheme with small quark masses. The corresponding 
scalar master integrals are evaluated using the methods 
and results of Refs. 
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461 ]. where different regulariza- 
tion schemes are translated into each other as described 
in Ref. 47[. Our treatment of rational terms of UV or 
IR origin is described in Appendix A of Ref. [H. More 
details on the two independent calculations of the virtual 
corrections in the gg channel will be presented elsewhere. 

In both evaluations of the real corrections the am- 
plitudes are calculated in the form of helicity matrix 
elements which have been generated with Madgraph 
4.1.33 48]. While the amplitudes for qq —> ttbbg have 
been checked with the spinor formalism of Ref. , those 
for gg — > ttbbg have been verified with an implementa- 
tion of off-shell recursion relations 5^, 51 , H3| ■ The sin- 
gularities for soft or collinear gluon emission are isolated 
via dipole subtraction 53, 54, 551 56] for NLO QCD cal- 
culations using the formulation [53j for massive quarks. 
One of the two calculations employs the automatic Mad- 
Dipole implementation of dipole subtraction [57j. After 
combining virtual and real corrections, singularities con- 



nected to collinear configurations in the final state cancel 
for "collinear-safe" observables after applying a jet algo- 
rithm. Singularities connected to collinear initial-state 
splittings are removed via MS QCD factorization by PDF 
redefinitions. In both evaluations the phase-space inte- 
gration is performed with multichannel Monte Carlo gen- 
erators [581 ] and adaptive weight optimization similar to 
the one implemented in Lusifer [59(. 

In the following we consider the process pp — > ttbb + X 
at the LHC, i.e. for i/s — 14TeV. For the top-quark 
mass, renormalized in the on-shell scheme, we take the 
numerical value m t = 172.6 GeV [H. All other QCD 
partons (including b quarks) are treated as massless par- 
ticles, and collinear final-state configurations, which give 
rise to singularities, are recombined into IR-safe jets 
using a fcT-algorithm |61|. Specifically, we adopt the 
fcT-algorithm of Ref. [62j and recombine all final-state b 
quarks and gluons with pseudorapidity |?7| < 5 into jets 
with separation y/ A<j> 2 + Ay 2 > D = 0.8 in the rapidity- 
azimuthal-angle plane. Requiring two b-quark jets, this 
also avoids collinear singularities resulting from the split- 
ting of gluons into (massless) b quarks. Motivated by the 
search for a ttH(H — > bb) signal at the LHC 0, we 
impose the following additional cuts on the transverse 
momenta, the rapidity, and the invariant mass of the 
two (recombined) b jets: px,b > 20 GeV, |j/b| < 2.5, and 
m bb > m b { 5CUt . We plot results either as a function of 
TO bb cut or f° r m b5 cut = 0- Note, however, that the jet 
algorithm and the requirement of having two b jets with 
finite PT,b in the final state sets an effective lower limit on 
the invariant mass m bb of roughly 20 GeV. The outgoing 
(anti)top quarks are neither affected by the jet algorithm 
nor by phase-space cuts. 



We consistently use the CTEQ6 [6J, l6J set of PDFs, 
i.e. we take CTEQ6L1 PDFs with a 1-loop running a s 
in LO and CTEQ6M PDFs with a 2-loop running a s in 
NLO, but the suppressed contributions from b quarks 
in the initial state have been neglected. The number 
of active flavours is Np — 5, and the respective QCD 
parameters are A^° = 165 MeV and A^ s = 226 MeV. In 
the renormalization of the strong coupling constant the 
top-quark loop in the gluon self-energy is subtracted at 
zero momentum. In this scheme, the running of a s is 
generated solely by the contributions of the light-quark 
and gluon loops. By default, we set the renormalization 
and factorization scales, /ir and /ip, to the common value 
/i = mt + m bbcut /2. 

In Figure [2] we show the scale dependence of the LO 
and NLO cross sections upon varying the renormaliza- 
tion and factorization scales in a uniform or an antipodal 
way. We observe an appreciable reduction of the scale 
uncertainty upon going from LO to NLO. Varying the 
scale up or down by a factor 2 changes the cross sec- 
tion by 70% in LO and by 34% in NLO. At the central 
scale, the full pp cross section receives a very large NLO 
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FIG. 2: Dependence of the LO and NLO cross sections of 
pp — > ttbb + X at the LHC for TO b g cut = and fj,o = rrh. 



correction of 77%, which is mainly due to the gluonic 
initial states. Introducing a veto on extra jets by requir- 
ing px, jet < 50 GeV reduces the iVfactor to roughly 1.2. 
For the qq channel we found a very small correction of 
2.5% [l|. The full LO and NLO cross sections are given 
by cj lo = 1488.8(1.2) fb and ctnlo = 2638(6) fb, where 
the numbers in parentheses are the errors of the Monte 
Carlo integration for 2 x 10 7 events. 

Figure[3]shows the LO and NLO cross sections as func- 
tion of the cut m bb cut on the bb invariant mass, where 
the bands indicate the effect from a uniform or antipo- 
dal rescaling of /xr. and //p by factors 1/2 and 2. The 
shown LO and NLO bands overlap in the whole consid- 
ered range in 77i bbcut , which is motivated by the search 
for a low-mass Higgs boson. In contrast to the pure qq 
channel [l[, the NLO corrections to the full pp process 
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FIG. 3: LO and NLO cross sections for pp — > ttbb + X at 
the LHC as function of m b g cut , with the bands indicating the 
scale dependence by varying fin, and /j,f by factors 1/2 and 2 
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J2). 



induce only a moderate distortion of the functional de- 
pendence on m bb cut . The reduction of the scale uncer- 
tainty from about ±70% to ±33% and the large impact 
of the NLO correction hold true for the considered range 
m m bb cut . 

In summary we have presented first results on the full 
NLO prediction for the process pp — > ttbb + X at the 
LHC. The NLO corrections appreciably reduce the un- 
physical scale dependence of the LO cross section, but at 
the same time enhance the cross section by a A"-factor of 
about 1.8 for the usual scale choice. This large correc- 
tion factor can be strongly reduced by imposing a veto 
on hard jets. It will be interesting to see how these NLO 
results influence the signal significance of tt-Higgs pro- 
duction with a Higgs boson decaying into a bb pair, to 
which direct ttbb production represents a serious back- 
ground. 

On the technical side the presented calculation con- 
stitutes the first complete NLO prediction for a hadronic 
process of the type 2 — > 4 particles. Speed and stability of 
the evaluation show good performance of the Feynman- 
diagrammatic approach that is augmented by dedicated 
reduction methods of tensor loop integrals for exceptional 
phase-space regions. 

On a single 3 GHz Intel Xeon processor, the evaluation 
of the virtual corrections for gg — > ttbb (including sums 
over colour and polarization states) takes about 160 ms 
per phase-space point. This remarkably high speed sug- 
gests that the employed reduction method might turn out 
to be a very efficient tool for various other multi-particle 
processes at the LHC. 
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